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We report the existence and properties of localized modes described by nonlinear Schrodinger 
equation with complex PT-symmetric Rosen-Morse potential well. Exact analytical expressions 
of the localized modes are found in both one dimensional and two-dimensional geometry with 
self-focusing and self-defocusing Kerr nonlinearity. Linear stability analysis reveals that these 
localized modes are unstable for all real values of the potential parameters although corresponding 
linear Schrodinger eigenvalue problem possesses unbroken T'T-symmetry. This result has been 
Mh. verified by the direct numerical simulation of the governing equation. The transverse power flow 

density associated with these localized modes has also been examined. 
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C^ . I. INTRODUCTION 

Recently, there has been considerable amount of attention to theoretical and experimental investigation 

CN , of light propagation in parity-time (VT) symmetric optical media Ill4l4l| . The interest in study of such 
> 
l/~j "PT-symmetric optical materials has its roots in quantum mechanics: the paraxial equation of diffraction is 

2_, • mathematically equivalent to that of quantum Schrodinger equation. Quantum mechanics requires that the 

^ ' spectrum of every physical observable should be real, which of course are satisfied by Hermitian operators. 

S ' However, Bender and Boettcher [l5|] pointed out that some non-Hermitian Hamiltonians with PT-symmetry 

cn can also exhibit an entirely real spectrum and may constitute unitary quantum systems without violating any of 

•• [ the axioms of quantum mechanics. Moreover, it has been shown that for a "PT-symmetric complex Hamiltonian, 

•'"! there may exist a threshold above which its eigenvalues are not real but become complex, and the system 

Vh . undergoes a phase transition because of spontaneous PT-symmetry breaking. In general the action of the 

- ■ ■ parity P and time T operators is defined as p — )• —p, x — )• —x and p — )• —p, x ^ x, i — > —i, respectively. A 

Hamiltonian with a complex PT-symmetric potential requires that the real part of the potential must be even 

function of position and the imaginary part should be odd [1G|]. In optics such complex PT-symmetric structure 

can be designed through a judicious designs that involve both optical gain/loss regions and the process of index 

guiding 4j-[6|]. In such settings, complex refractive index distribution n{x) = nii{x) + i nj{x) plays the role of an 

optical potential so that the index guiding nR{x) and the gain/loss profile ni{x) satisfy nji{x) = nji{—x) and 

ni{x) = —ni{—x), respectively. Unusual exotic phenomena like VT phase transition, band merging, double 

refraction, non-reciprocity [a, |20|, |2l|] , and unidirectional invisibility 17H19l| etc have been reported to exist in 



'Electronic address: |bikash.midya@gmail.com| 
^Electronic address: |rroychoudhuryl23@gmail.com| 



linear "PT-symmetric complex optical media. Spontaneous "PT-symmetry breaking has been experimentally 
observed in active or passive VT dimers [3, \^ and periodic lattices [l|]. These findings, in turn, have stimulated 
considerable research activity in the non-linear PT-symmetric systems as well. 

In the nonlinear domain, a novel class of one and two dimensional localized modes were found to exist 



below and above the phase transition point 22] and the interplay between the Kerr nonlinearity and the VT 



threshold was investigated |6|]. Subsequently, nonlinear modes are studied in complex PT-symmetric periodic 
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32l |. Gaussian [33;], Bessel |23|], Scarf-ll |3l|], and harmonic [35|| potentials, as well as in a harmonic trap 
with a rapidly decaying PT-symmetric imaginary component [34|]. Stable localized modes in a PT-symmetric 



slab waveguide with distributed gain and loss are found in 36|]. Existence of optical solitons in PT-symmetric 



nonlinear couplers with gain/loss [371,1381], gap solitons in PT-symmetric optical lattices j32] and optical defect 
modes in PT-symmetric potentials [40] are also reported. Stable ID and 2D bright spatial solitons are found 



to exist in defocusing Kerr media with PT-symmetric Scarf II potentials [24|]. Also, it has been found that 
the gray solitons in PT-symmetric potentials can be stable [41]. However the existence of nonlinear localized 
modes in yet another important potential e.g. complex PT-symmetric Rosen-Morse well has not been reported 
so far. The complex Rosen-Morse potential well is characterized by the same real component as the complex 
Scarf-ll potential, however, its imaginary component is different. In fact, in contrast with the real component, 
the imaginary potential component doesn't vanish asymptotically, rather it tends to a finite value. This is the 



reason why the phenomenon of spontaneous breakdown of PT-symmetry is elusive in such system [25|, l26l ]. 
Nevertheless, the bound state energy eigenvalues of PT-symmetric Rosen-Morse potential well undergoes a shift 
from negative to positive domain for certain range of parameters which controls the strength of the potential. 
In this paper, we investigate the propagation of nonlinear beam in a single VT waveguide cell which 
is characterized by the nonlinear Schrodinger equation with complex Rosen-Morse potential well. Specifically, 
the existence of the spatial localized modes have been reported in both one-dimensional and two-dimensional 
settings with self-focusing and self-defocusing Kerr nonlinearity. We have shown, with the help of linear stability 
analysis of the one-dimensional localized modes, that though the spontaneous breakdown of PT-symmetry does 
not occur in complex Rosen-Morse well the localized modes corresponding to nonlinear Schrodinger equation 
are always unstable. This linear instability has been verified by direct numerical simulation of the governing 
equation. The transverse power flow density associated with these nonlinear localized modes has also been 
examined. 

II. LOCALIZED MODES IN PT-SYMMETRIC COMPLEX ROSEN-MORSE WELL 

A. Mathematical Model 

We consider optical wave propagation in a Kerr nonlinear PT-symmetric potential. In this case, (1 -|- 
l)-dimensional optical beam propagation along longitudinal z direction is governed by the following non-linear 
Schrodinger like equation [a lid] 

i-Q^ + -Q^ + [^(^) + ^^(2^)] ^ + ^1*1 * = 0- (2.1) 



Here ^(x,z) is slowly varying complex electric field envelop, x is the transverse co-ordinate, and u = ±1 
represent the self-focusing and self-defocusing nonlinearity respectively. V{x) and W{x) are the real and 
imaginary parts of the complex PT-sym metric potential such that V{—x) = V{x) and W{—x) = —W{x). 
Physically, V{x) is responsible for the bending and slowing down of light, and W{x) can lead to either 
amplification (gain) or absorption (loss) of light within an optical material. 

The optical beam propagation in a single VT cell is important to understand light self-trapping in complex 
optical lattices. In order to investigate the optical beam propagation in a single VT cell, we consider the 
complex ■PT-symmetric Rosen-Morse potential well as 

V{x) = -a{a + l)sec\)?x, 
W{x) = 2ht&nh.x, 

where a and h characterizes the strength of the real and imaginary parts of the potential, respectively. Both the 
real and imaginary part of this potential are shown in figure[T^ for the potential parameters a = .75 and h = .8. 



The linear Schrodinger eigenvalue problem for the potential 12.21 has been thoroughly studied in [25|, 129]. It has 
been shown that all the bound state energy eigenvalues corresponding to the linear Schrodinger equation of 
complex Rosen-Morse well are real so that the spontaneous breakdown of 'PT-symmetry never occur. However, 
the energy eigenvalues 

52 

An = -(a - nf + -^, n = 0, 1, ...rimax < a. (2.3) 

[a — nj^ 

begin to shift from negative to the positive domain when the strength of the non-Hermiticity is increased. In 
fact all the energy eigenvalues become positive whenever ^y\b\ > a^. Here, we search for stationary solution of 
the nonlinear equation ()2.ip in the form ^(x,z) = (j){x)e'^^'^ , where A is the real propagation constant, and the 
complex function (l){x) satisfies following equation 

—I - [a(a -h 1) sech^ X - 2ibt<mhx\ (j) + a\(t)\^ (j) = X(j). (2.4) 

In the following we report the existence and linear stability of the localized modes of the above nonlinear equation 
(|2.4p for both the self- focusing and self-defocusing cases. 

B. Analytical solutions and their linear stability 

1. Self-focusing case (a = 1) 
For a = 1, equation (j2.4p is found to admit an exact analytical expression of the localized mode of the form 



(p{x) = \/a2 + a + 2 sechx e''"', (2.5) 

where A = 1 — 6^. In figure [1^, the real and imaginary parts of this spatial soliton have been shown for a = .75 
and 6 = .8 and A = .36. To focus on the properties of this non-linear solution, we examine following three 
quantities: the transverse power flow density (Poynting vector) S across the beam, the power P, and the 
linear stability of these localized modes. For the nonlinear modes given in equation (12. 5p . the Poynting vector 
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FIG. 1: (Color online) (a) Real and Imaginary parts of the Rosen-Morse Potential Well; (b) Real and Imaginary 
parts of the localized modes (t>{x) in the self-focusing medium; (c) The transverse power (Poynting vector) S{x); (d) 
Power P{a) as a function of potential parameter a. In (a), (b) and (c) we have considered a = .75, 6 = .8, A = .36 
and (7 = 1. 



S = ^{(t)4>x ~ 4'*4>x) = b{a'^ + a + 2) sech^ X depends on the sign of the strength, b, of the imaginary part of 
the potential. It may be both negative and positive for negative and positive values of b respectively. However, 
we consider only positive values of b in which case S is everywhere positive and the power flow in the VT cell 
is in one direction, i.e. from the gain towards loss domain. For a = .75, b = .8, the transverse power flow is 
shown in figure [Tfc). For the localized modes (j2.5p the power P is calculated as 



P{a,b) 



{x)\'^dx = 2{a'^ + a + 2). 



(2.6) 



Clearly, the power is independent of the parameter b. In figure [T{d) we have shown the power P as a function 
of the real potential strength parameter a. It remains positive for all values of the parameter a. The Power 
increases with the increasing absolute value of the parameter a and becomes minimum for a = —.5. At a = —.5, 
the amplitude of the localized modes is also minimum. 

In order to determine the linear stability properties of the self-trapped localized modes obtained here, we 
consider small perturbation to the solution "i>{x,z), of the form 
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^{x, z) = 0(x)e*^^ + [[f{x) + g{x)] e"^ + [/*(x) - g*{x)] e"*^} 



A\z 



where f{x) and g{x) are infinitesimal perturbation eigen-functions such that |/|, \g\ <^ 
perturbation growth rate. By linearizing the equation (j2.ip around the localized solution (/ 
functions / and g satisfy the following eigenvalue problem 






(2.7) 

■q stands for the 
, we find that the 



(2.8) 



where ti = d^x + {V + iW) + a\4>\'^ - A and £2 = d^x + {V + iW) + 3cr|(/>p - A. The linear stability of 
the localized modes (/>(x) depends on the nature of the eigenvalue 77. The "PT-symmetric nonlinear localized 
mode is unstable if r/ has any positive real part, because for 3?(r/) > the corresponding perturbed nonlinear 
eigenmodes (12. 7p would grow exponentially with z. The eigenvalues rj can be obtained by solving equation (12. 8p 
with the help of several numerically techniques [29]. In this paper we have used Fourier collocation method 



301 1 . Our numerical investigations corresponding to the nature of the eigenvalue r] reveal that i] never becomes 
purely imaginary. It has always non vanishing positive real part for all real values of potential parameters a 
and b. This implies that the nonlinear localized modes obtained here are always unstable. The results of linear 
stability analysis are corroborated by direct numerical simulations of Eq. (j2.ip using the solution (j2.5p as initial 
condition i.e. "^{x,0) = (j){x). In figure [2] we have shown the localized modes cl){x), unstable intensity profiles 
|^(x,2;)p and corresponding linear stability spectra for very small and large values of the parameter b. As 
expected, localized modes which are predicted to be unstable fail to maintain their original shapes. The reason 
behind such instability of the localized modes is that unlike the real part, the imaginary part of the Rosen-Morse 
potential well does not vanish asymptotically. Therefore gain/loss remains in the system even far from the place 
of localization and any small fluctuations of the field is amplified/absorbed, eventually leading to instability. 
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FIG. 2: (Color online) (a) Plots of the real (solid blue curve) and imaginary (dotted-dashed red curve) parts of 



{x) and 



(solid black curve) in self focusing medium; (b) The evolution of field intensity |^(a;,z)p (c) 



Numerically computed stability spectra corresponding to the figure [2^. In all these cases we have considered a = 
.1, b = .03, A = .999 and cr = 1. (d), (e) and (f) Plots of the same quantities as in figure[2{a), (b) and (c) respectively, 
for a = .1, 6 = 3, A = —8, a — I. 



It is worth mentioning here that unlike the "PT-symmetric Scarf II potential (for which the solitons are 
stable below the certain critical value of the imaginary potential component and become unstable above this 
critical value [a, |3lt|), localized modes in the complex Rosen-Morse potential well, discussed here, are linearly 
unstable for all real values of non-Hermiticity parameter b. This result is valid in spite of the fact that the 



Hamiltonian corresponding to the linearized version of equation (J2.4I) possesses unbroken "PT-symmetry (all 

energy eigenvalues are real). Nevertheless, the parameters range y^ > a^, for which the transition from the 

negative to positive energy corresponding to the linear Schrodinger eigenvalue problem takes place, does not 
affect the instability of the localized modes. Only the magnitude of the localized modes differs. 



2. Self-defocusing case (a = —1) 



Optical beam propagation in nonlinear self-defocusing Kerr medium is governed by the equation (j2.ip for 
(7 = — 1 and its corresponding stationary solutions satisfy equation (j2.4p . Like the self focusing case, here 



equation (12. 4p admits an exact solution (j){x) = a/— (c?-|-a + 2) sechx e*^^. Note that these non-linear modes 
in self-defocusing case are very similar to those in the self-focusing case. Only the amplitude of the localized 
modes are different. Nevertheless, like the self-focusing case here also the localized modes are unstable for all 
o, 6 E K. In figures [SJa), and[3Jb), we have shown such nonlinear modes and corresponding unstable intensity 
evolution for the parameter values a = 1, 6 = .4, A = .84 and a = —1. Numerical solution of the eigenvalue 
problem (|2.8p has been plotted in figure [SJc) which also implies that the the corresponding modes are linearly 
unstable. 
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FIG. 3: (Color online) (a) Plots of the real (solid blue curve) and imaginary (dotted-dashed red curve) parts of 
the localized modes 4>{x) in self-defocusing media; (b) Unstable intensity evolution |^(x, z)p corresponding to the 
figure [3^; (c) Plot of the corresponding stability spectra obtained numerically. In all these cases we have considered 
potential parameters a = 1, 6 = .4, A = .84 and cr = — 1. 



III. LOCALIZED MODES IIS! TWO-DIMENSIONS 



Finally, we discuss the formation of nonlinear localized modes in two-dimensional Rosen-Morse potentia 
The two-dimensional generalization of the equation ()2.ip . with the self-focusing nonlinearity, is given by 22] 



Z----F V^^-F [V{x,y) + iW{x,y)]^! + |^p^ = 0, 



(3.1) 



where V 



dx' 



+ 



dip 



is the two-dimensional Laplacian. The two-dimensional complex Rosen-Morse potential 



well, which obey the "PT-symmetric requirements V{—x,—y) = V{x,y) and W{—x,—y) = —W{x,y), can be 
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FIG. 4: (Color online) (a) Real part of the 2D Rosen-Morse potential, (b) Imaginary part of the 2D Rosen-Morse 
Potential; (c) Plot of the 2D localized modes |(?!)(a;,y)p. (d) Transverse power flow vector S indicating the power flow 
from gain towards loss regions. In all these cases we have considered a = 1.25 and 6 = .5, A = 1 and a = \. 



considered as 

y(x, y) = 2(sech^ x -|- sech^ y) — [c? -|- a -|- 2) sech^ x sech^ y 
W{x^y) = 46(tanha; -|- tanhy). 

The stationary solutions of the equation (13. ip can be assumed in the form 

^{x,y,z)=4>{x,y)e'^'+''^^'y^ 

where <j){x,y) and the phase 9{x,y) satisfy the following two equations 

V^(j) - \Ve\^(j) + V{x, y)cp + </.3 = A(/), 
(t)W^e + 2V0. V(/. + W{x, y)(p = 0. 

respectively. 

A nonlinear solution to equation p.4p that satisfies (/)—)• as {x,y) — > ±oo is obtained as 



(3.2) 



(3.3) 



(3.4) 



(l){x, y) = V a^ + a + 2 sech x sech y, 



(3.5) 



with the phase 6{x,y) = b{x + y) and the propagation constant is given by A = 2 — 46^. Figures |4{a),(b) 
show the real and imaginary parts of the 2D Rosen-Morse potential well. Two-dimensional soliton \(l){x,y)\^ is 
shown in figure HI^c). In all these cases we have considered a = 1.75 and a = 1. To understand the internal 
structure of the two dimensional self-trapped modes, we calculate the two-dimensional transverse power flow 
vector S = b{a? + a + 2) (sech^ x,sech^ y). In figure |4{d), we have shown such 2D transverse power flow for 
a = 1.25,6 = .5A = 1,0" = 1, which indicates energy exchange from gain towards loss regions. 

IV. SUMMARY 



To summarize, we have investigated the existence and properties of nonlinear localized modes in a single 
VT waveguide cell characterized by the nonlinear Schrodinger equation with complex Rosen-Morse potential 



well. The closed form expressions for the localized modes in such one- and two-dimensional self-focusing and 
self-defocusing Kerr nonlinear media are obtained. The transverse power flow density is shown to remain 
positive for some parameter values which indicates that power flow is in a single direction, mainly from gain 
towards loss regions. However, linear stability analysis of the one-dimensional solitons reveals that these solitons 
are unstable over the whole range of the potential parameter in spite of the fact that corresponding linear 
Schrodinger eigenvalue problem possesses unbroken PT-symmetry. The main reason behind such instability 
is that unlike the real part, the imaginary part of the complex Rosen-Morse potential well does not vanish 
asymptotically. Therefore any small fluctuation in the field intensity is amplified (or absorbed) which leads to 
the instability. The results presented here definitely encourage one to search for the stable localized modes (if 
any) of the nonlinear Schrodinger equation with the "PT-symmetric Rosen-Morse potential well in the presence 
of higher-order/competing or other nonlinearities. 
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